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The systematic investigation of hyperdeformation (HD) at high spin in the Z = 40 — 58 part of 
the nuclear chart has been performed in the framework of the cranked relativistic mean field theory. 
The properties of the moments of inertia of the HD bands, the role of the single-particle and necking 
degrees of freedom at HD, the spins at which the HD bands become yrast, the possibility to observe 
discrete HD bands etc. are discussed in detail. 



o 
o 

(N 



PACS numbers: 21.60.Jz, 27.50.+e, 27.60.+j, 21. 10. Ft, 21.10.Ma 



> 

0\ 

o 
o 

o 

O 



X 



I. INTRODUCTION 

Since the discovery of superdeformation (SD) in -'^^^Dy 
two decades ago nuclear SD has been in the focus 
of attention of the nuclear structure community; it has 
been discovered in different mass regions and extensively 
studied experimentally ^ and theoretically (see, for ex- 
ample, Refs. [1,13,0] and references therein). New phe- 
nomena such as identical bands [sl] were discovered, and 
rich variety of experimental data allowed to test mod- 
ern theoretical tools under extreme conditions of large 
deformation and fast rotation. 

It was known for a long time from harmonic oscilla- 
tor studies that even more elongated shapes, called 
as hyperdeformed (HD) and characterized by the semi- 
axis ratio of around 3:1, are possible. The existence of 
such stable shapes was later confirmed in the macro- 
scopic+microscopic (MM) method 0, H, H, 0, HH, [H, 
[Tsl. IT^. [Tst. Theoretical results on the states located in 
third (HD) minima are also available in self-consistent 
Hartree-Fock-|-Bogoliubov (HFB) approaches based on 
the Skyrme and Gogny forces (see Refs. [H, [O, [ill and 
references quoted therein), and relativistic mean field ap- 
proach [19]. However, these results are restricted to spin 
zero states, which are difficult to measure in experiment. 
To our knowledge, the description of the HD states at 
high spin within the self-consistent approach has been at- 
tempted only in ^°*Cd 20] [within the cranked relativis- 
tic mean field (CRMF) method] and in four ^ ~ 40 mass 
nuclei [2l| [within the cranked Skyrme-Hartree-Fock ap- 
proach]. The general feature of all these calculations is 
the fact that the semi-axis ratio of the HD shapes is less 
than 3:1 [1]. 

Let us mention two examples of such studies: one at 
spin zero, another at high spin. In actinide nuclei, the 
HD states are so-called third minima states around ^'^^Th 
[ill . [23l . [2^. In these nuclei, the second saddle point is 
split, leading to the excited reflection-symmetric and re- 
flection asymmetric configurations with large quadrupole 
and octupole deformations, /32 ^ 0.9 and (3^ ~ 0.35. The 
density distribution at the HD minimum resembles a di- 
nucleus consisting of a nearly-spherical nucleus around 
the doubly-magic nucleus -'^^^Sn and a well-deformed frag- 



ment from the neutron-rich A ^ 100 region [ll|. Unfor- 
tunately, it is very difficult to study the HD states at 
low spin in experiment. In order to overcome this prob- 
lem, one should use the fact that the larger moment of 
inertia connected with the larger deformation drives the 
nucleus towards larger deformations with increasing an- 
gular momentum; the HD minimum is thus favored by ro- 
tation and becomes ultimately yrast at high spin. For ex- 
ample, cranked Nilsson-Strutinsky calculations suggested 
the existence of very elongated high-spin minima in nu- 
clei around ^^^Yh ^. These HD bands are expected to 
become yrast at spin around SOft. 

On the experimental side, very little was known about 
hyperdeformation apart from some indications of this 
phenomenon at low spin in the uranium nuclei [25| and 
light nuclei like ^^C [2g] and the observation of the HD 
ridge structures at high spin in the A ~ 150 mass region 
[271 . i28i] . Recent observation of the very extended shapes 
in ^°*Cd [29l . [30l . strongly motivated by earlier calcula- 
tions of Ref . and more recent studies of Ref . ^15*] , has 
renewed interest in the study of hyperdeformation at high 
spin. Although the hyperdeformed nature of the bands 
in this nucleus has not been confirmed in the subsequent 
cranked relativistic mean field analysis of Ref. (see 
also Sect. VB in Ref. [3l|), this experiment provided a 
strong motivation for subsequent experimental searches 
in the A ~ 125 mass region (see Refs. [13, [13, [11]) and 
theoretical studies of Refs. 'B', '3T] within the framework 
of the MM method. These experiments revealed rota- 
tional patterns in the form of ridge-structures in three- 
dimensional (3D) rotational mapped spectra with dy- 
namic moments of inertia J*-^^ ranging from 63 to 111 
MeV^^ in 12 different nuclei [l^]; the values around 110 
MeV-i observed in "^Xe, i24xe and 124,125^3 suggest 
that the HD structures were populated in these experi- 
ments. However, no discrete rotational HD bands have 
been identified. It is also necessary to mention that sev- 
eral previous atternpts to search for high spin HD struc- 
tures in "7Gd [11[3|, i52Dy [13, [13, and i^syb 36] did 
not lead to convincing evidences for discrete HD bands. 

So far, theoretical investigations of HD at high spin 
were carried out mainly in the framework of the MM 
method. One of the main goals of the current manuscript 
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FIG. 1: (Color online) The chart of nuclei in the Z = 40 — 58 region. Only experimentally known nuclei are shown. Experimental 
data on superdeformed nuclei are taken from Ref. 0|. The nuclei in which the search for HD structures has been performed in 
the HLHD experiment are taken from Ref. [2^ . 



is to perform for the first time a systematic study of HD 
within the framework of fully self-consistent theory, the 
CRMF theory. Fig. [1] shows the part of the nuclear chart 
where our studies are performed. We restrict our investi- 
gation to even-even nuclei; the only exceptions are odd- 
mass nuclei ^^^I [in which extremely SD doubly magic 
band has been found] and ^^^'-^^"^Xe, ^^'^I and ^^^Cs [which 
are used in the study of the relative properties of the HD 
bands] . In each isotope chain we consider nuclei ranging 
from the most proton-rich ones up to the ones located at 
the neutron-rich side of the /3-stability valley. Neutron- 
rich nuclei beyond the valley of the /3-stability are ex- 
cluded from consideration because of the experimental 
difhculties of studying them at high spins relevant for 
HD. With the goal to guide future experimental explo- 
rations and to find the nuclei in which the HD may be 
studied with current and future experimental facilities, 
we define the spins at which the HD bands become yrast 
in these nuclei. In addition, available experimental data 
on the HD ridge-structures in the Te, Xe, and Cs nuclei 
are analyzed. The general features of the HD bands are 
outlined. 

The role of the single-particle degrees of freedom at 
hyperdeformation has not been studied in detail till now. 
One of the major goals of the current manuscript is the 
study of their role, and it is motivated by the desire to un- 
derstand to what extent theoretical methods developed 
in the study of the SD bands are also applicable to the 
HD bands. It is very unlikely that the spins, parities and 
excitations energies of the HD bands will be known in 
the initial stage of their experimental study. The direct 
test of the structure of the wave functions of the single- 
nucleonic orbitals (e.g. via magnetic moments) will also 
not be possible at that sta ge. Thus, similar to the case 
of superdeformation d, [btI. Issl. [39| , the relative proper- 
ties of different HD bands may play an important role in 
the interpretation of their structure. In this context, it 



is important to understand which changes of the single- 
particle orbitals are involved in going from one HD band 
to another, and how they affect physical observables like 
dynamic moments of inertia J^^' , transition quadrupole 
moments Qt, total spin /, etc. In particular, we will 
study whether the theoretical methods which were sys- 
tematically used in the configuration assignment of the 
SD bands are also applicable to the HD bands. These in- 
clude the methods based on the relative properties of the 
dynamic moments of inertia J*-^^ [1, , on the effective 
alignments ief / 0, Issi . js^ and on the relative transition 
quadrupole moments AQt [iO, Slj . 

The manuscript is organized as follows. The definition 
of physical observables and the details of numerical cal- 
culations are discussed in Sect. |TT1 The spins at which 
the HD bands become yrast, the regions of nuclear chart 
where the experimental search for the HD structures may 
be successful and the general properties of the HD bands 
are outlined in Sect. lIIII The data obtained in the search 
of the HD structures in the yl ^ 120 mass region and the 
single-particle degrees of freedom are also analysed in this 
section. Sec. IIVI is devoted to the analysis of extremely 
superdeformed (ESD) structure in ^^^I. The calculations 
predict the existence of doubly magic ESD structure in 
this nucleus with the deformations being close to HD, 
which may be observed with the current generation of 
7-ray detectors. Finally, Sect. fVl contains the main con- 
clusions of our work. 



II. THE DETAILS OF THE CALCULATIONS 

In the relativistic mean field (RMF) theory the nu- 
cleus is described as a system of point-like nucleons, 
Dirac spinors, coupled to mesons and to the photons 
[13, m, 113]. The nucleons interact by the exchange of 
several mesons, namely a scalar meson a and three vec- 
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tor particles cu, p and the photon. The cranked relativis- 
tic mean field (CRMF) theory 0, IH, [H, [13 represents 
the extension of RMF theory to the rotating frame. It 
has successfully been tested in a systematic way on the 
properties of different types of rotational bands in the 
regime of weak pairing such as normal-deformed [48l | , su- 
perdeformed P.|49|| as well as smooth terminating bands 

In the current study, we restrict ourselves to reflec- 
tion symmetric shapes since previous calculations in the 
MM method show no indications that odd-multipole (oc- 
tupole, ...) deformations play a role in the SD and HD 
bands of the nuclei covered by our study and in the 
HD bands of the A 110 — 125 [s^l mass region. 



A. Physical observables 

Similar to the case of the SD bands, it is reasonable 
to expect that the HD bands will not be linked to the 
low-spin level scheme for a long period of time. Thus, the 
spins and parities of the HD bands will not be known and 
it will not be possible to define the kinematic moment of 
inertia J^^^ since it depends on the absolute values of the 
spin. In such a situation, the dynamic moment of inertia 
J*^^^ will play an important role in our understanding of 
the structure of the HD bands. This is similar to the case 
of the SD bands (see Refs. d, [13] )■ Other observables, 
such as transition quadrupole moments Qt ^-nd effective 
(relative) alignments ieff, will also be important. 

In the CRMF calculations, the rotational frequency 
Qx, the kinematic moment of inertia J^^^ and the dy- 
namic moment of inertia J^^-* are defined by 
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The charge quadrupole Qo a-nd mass hexadecupole Q40 
moments are calculated by using the expressions 
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go = eJ—^{r^Y^,)l + 2{r^Y^^)l, (4) 



moments using simple relation [1^ |53, [5 
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where R — \.2A^/'^ fm is the radius of the nucleus, and 
is a quadrupole moment of the X-th (sub)system ex- 
pressed in fm^. Here X refers cither to proton {X = Z) 
or neutron {X — N) subsystem or represents total nu- 
clear system {X = A). This expression, however, neglects 
the higher powers of fj^ and higher multipolarity defor- 
mations /?4,/36,--- [131, which play an important role at 
hyperdeformation. 

Considering that the definition of the deformation is 
model dependent [s^ . and that this quantity is not 
experimentally measurable, we prefer to use transition 
quadrupole moment Qt for the description of deforma- 
tion properties of hyperdeformed states. This is experi- 
mentally measurable quantity, so in the future our pre- 
dictions can be directly compared with experiment. The 
deformation properties of the yrast SD band in ^^^Dy 
(which is one of the most deformed SD bands [1^) are 
used as a reference. This is done by introducing normal- 
ized transition quadrupole moment Qi°™{Z^A) in the 
{Z, A) system 
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Q40 = (r4r4o)„ + (r"no)„ 



(5) 



This equation is based on the ratio 

gnorm(^^^)/g^(152j3y) calculated USiug Eq. ® 

under the assumption that the /32-values in the (Z, A) 
system and in ^^^Dy are the same. We use the value 
(5t(^^^Dy) = 18.73 eb obtained in the CRMF cal- 
culations with the NLl parametrization of the RMF 
Lagrangian for the yrast SD band in ^^^Dy at / = 60?i 
in Ref. [J]. Thus, in first approximation (neglecting 
the higher powers of /32 and higher multipolarity de- 
formations PijPe,...)) the equilibrium deformation of 
the band in the {Z,A) system having the QJ'°™(Z,yl) 
value is the same as in the yrast SD band of ^^^Dy. 
We describe the band as hyperdeformed if its Qt value 
exceeds QJ'°™(Z,^) by at least 40%. This criteria 
is somewhat relaxed in the Z — 40, 42, 44 nuclei for 
which the band is defined as HD if its Qt value exceeds 
Q^'"-"'{Z,A) by at least 30%. 

The effective (relative) alignment «e// between two 
bands is defined as the difference between the spins of 
two levels in bands A and B at the same rotational fre- 
quency rix (ssj : 



where the labels p and n are used for protons and neu- 
trons, respectively, and e is the electrical charge. At ax- 
ially symmetric shapes, typical for the hyperdeformed 
states, the transition quadrupole moment Qt is equal to 
Qo- 

The quadrupole deformation for axially-symmetric 
shapes is frequently defined in self-consistent calcula- 
tions from calculated and/or experimental quadrupole 
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This quantity has been used frequently in the analysis 
of the single-particle structure of the SD bands and the 
configuration assignment (see Refs. [38l [39| and refer- 
ences quoted therein) . It depends on both the alignment 
properties of the single-particle orbitals(s) by which the 
two bands differ and the polarization effects induced by 
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the particles in these orbitals [5l| 
related to nuclear magnetism. 



The latter are in part 



Because the pairing correlations are relatively weak in 
the HD bands of interest (see Sect. IIII C|) . their intrin- 
sic structure can be described by means of the dominant 
single-particle components of the hyperintruder states oc- 
cupied. The calculated configurations will be labeled by 
[p,nin2\, where p, ni and n2 are the number of proton 
N = 7 and neutron N = 7 and N — 8 hyperintruder 
orbitals occupied, respectively. For most of the HD con- 
figurations, neutron N = 8 orbitals arc not occupied, 
so the label n2 will be omitted in the labeling of such 
configurations. 

Single-particle orbitals are labeled by [A^n^Ajri**^". 
[NnzA]n are the asymptotic quantum numbers (Nilsson 
quantum numbers) of the dominant component of the 
wave function at il^ — 0.0 MeV. The superscripts sign 
to the orbital labels are used to indicate the sign of the 
signature r for that orbital (r = ±i). 

The spins at which the SD and HD configurations be- 
come yrast in the calculations are defined as crossing 
spins and , respectively. 
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FIG. 2: (Color online) Potential energy surfaces (PES) ob- 
tained in the axially symmetric RMF calculations without 
pairing in the ^''^Ce nucleus. The results of calculations with 
Nf = 12, 14 and 26 are shown. In all these calculations, Nb 
is fixed at 26. The results with Nf = 26 correspond to a fully 
converged solution: the binding energies do not change with 
further increase of Nf- The gaps in the PES lines are due to 
jumps of the solution from one single-particle configuration 
to another. The same single-particle configurations are com- 
pared at the same value of charge quadrupole moment. The 
normalized value of transition quadrupole moment Q""'''" cor- 
responding to the deformation of the yrast SD band in ^^^Dy 
is indicated by arrow. The range of hyperdeformation is also 
indicated. 



B. Numerical scheme of the CRMF calculations 

The CRMF equations are solved in the basis of 
an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates characterized by the deformation 
parameters Po and 7 and oscillator frequency Hloq = 
MeV, for details see Refs. 0, |4|. The trun- 
cation of basis is performed in such a way that all states 
belonging to the shells up to fermionic Np and bosonic 
Nb are taken into account. 

The impact of the truncation of basis on the numer- 
ical accuracy of the calculations has first been studied 
in the axially symmetric RMF code, see Fig. [2l In the 
mass region of interest, the calculations with Np — 12 
provide a reasonable approximation to the fully conver- 
gent Np = 26 solution up to a deformation typical for 
the SD shapes. However, this truncation scheme be- 
comes a poor approximation when the quadrupole mo- 
ment appreciably exceeds the one corresponding to the 
lower limit of HD; the difference between the Np = 12 
and Np = 26 solutions increases rapidly with the in- 
crease of quadrupole moment (see Fig. [5]). On the other 
hand, in this quadrupole moment range the results of the 
calculations with Np = 14 are closer to exact solution, 
although still exceeding it by ~ 1 — 2 MeV at the up- 
per end of the calculated quadrupole moment range. It 
was tested that with the decrease of the mass, the differ- 
ence between the Np = 14 and Np = 26 solutions will 
also decrease as well, so that the difference falls within 
the range of 1 MeV for the majority of the nuclei under 
study. 

These conclusions have also been tested in triaxial 
CRMF calculations. It was concluded that physical ob- 
servables of interest are described with sufficient numer- 
ical accuracy when Np = 12 is used for the SD and ND 
states and Np = 14 for the HD states. Thus, we em- 
ploy a hybrid calculational scheme in which the CRMF 
solutions in the ND- and SD minima are sought us- 
ing Np = 12, while the ones in the HD minima using 
Np = 14. In aU CRMF calculations, we use Nb = 20. 
In order to eliminate the numerical inaccuracies in the 
definition of the crossing spin , the yrast ND/SD 
configurations, which are crossed by the yrast HD con- 
figuration, were recalculated in the crossing region using 
Np — 14, and only then the crossing spin was defined. 
One should keep in mind that even with Np = 14 the 
spins at which the HD configurations become yrast in 
the calculations may be overestimated by 1 — 2ft when 
the deformation of the HD configurations exceeds appre- 
ciably the one corresponding to the lower limit of HD. 

When searching for different types of rotational struc- 
tures it is important to find the solutions in all local min- 
ima which are close to the yrast line in order to properly 
define the crossing spins between the rotational struc- 
tures of different nature. This is easily achievable in 
the macroscopic-|-microscopic approach by creating po- 
tential energy surfaces (PES) in the deformation s pac e 
covering quadrupole and triaxial deformations |55||. 
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FIG. 3: (Color online) The same as in Fig. [2l but for the 
results obtained in the axially symmetric RMF calculations 
with pairing using different parametrizations of the RMF La- 
grangian and Nf = 26. Figure shows the binding energies 
normalized with respect to the lowest energy of the lowest 
potential energy curve. 
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FIG. 4: (Color online) The crossing spins (left panels) at 
which the SD (solid circles) and HD (open squares) configura- 
tions become yrast and their transition quadrupole moments 
Qt (right panels) for the Ce, Ba and Xe isotopes. The values 
for the SD configurations are shown only in the cases when 
they become yrast at lower spins than the HD configurations. 
The normalized transition quadrupole moments Q""'""' corre- 
sponding to the deformation of the yrast SD band in ^^^Dy 
are also shown. 



However, the computational cost to create similar PES 
in the self-consistent models is enormous, thus, it has 
never been attempted in rotating nuclei. In order to 
overcome this problem, we use the fact that in self- 
consistent approaches without pairing the deformation 
of the basis defines to a large extent the local minima 
where the solutions will be obtained. Thus, the so- 



lutions in the ND minima, including triaxial ones, are 
searched using three combinations of the deformation of 
basis: (/3o = 0.30,7 = -30°), (/3o = 0.30,7 = 0°), and 
{/3q — 0.30,7 — +30°). In a similar way, the solutions in 
the SD minima are searched using the following combina- 
tions of the deformations of basis (/3o = 0.65, 7 = —30°), 
(/3o = 0.65,7 = 0°), {Po = 0.65,7 = +30°), and 
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(/3o = 0.8,7 = 0°)- The latter deformation of basis 
also leads frequently to the HD solutions. The defor- 
mation of basis {Po = 1.0, 7 = 0°) has been used for the 
search of the solutions in the HD minima. Non-zero 7- 
deformations of basis at large /3o lead either to the same 
solution as 7 = 0° or to the highly excited configurations. 
For each of the above mentioned values of the deforma- 
tion of basis, the lowest in energy solutions are calculated 
as a function of spin, and the yrast line is formed from 
these solutions. 
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regime in the SD and ND minima 0, [H, IH, Ull, the 
single-particle energies for the nuclei around the valley 
of f3 stability [39l. l53 | and the excitation energies of the 
SD minima j57l|. NL3 [SS^ is an alternative parametriza- 
tion, the quality of which has been tested in rotating 
nuclei (but less extensively than in the case of NLl) 
[sol , [isl . I49L [59I I . Some results with this parametrization 
will be presented. Few results obtained with the NLSH 
60] and NLZ [6l| parametrizations will be shown in Sect. 
IIII CI in order to illustrate the possible spread of calcu- 
lated quantities. It is necessary to keep in mind that the 
quality of the NLSH parametrization in respect of the 
description of rotational properties of the nuclei as well 
as their single-particle energies is not as good as that of 
the NLl and NL3 ^Si,!!!], and the force NLZ has not 
been tested in that respect. 

The spins at which the rotational structures belong- 
ing to different minima in potential energy surfaces be- 
come yrast depend in general on the relative energies 
of these minima and on the moments of inertia of rota- 
tional structures in these minima. Previous experience 
shows that different parametrizations of the RMF La- 
grangian give similar moments of inertia for the same 
configuration [H, 0, [H, [111 (see also Fig. [H below). 
Fig. [3] also illustrates that the potential energy surfaces 
at spin zero as a function of charge quadrupole moment 
obtained with the NLl and NL3 parametrizations are 
similar in shape. These two facts suggest that the HD 
configurations should become yrast at approximately the 
same spins in both parametrizations: this conclusion is 
confirmed in Sect. IIII Al It is interesting to note that the 
NL3 curve in Fig. [3] is similar to the one obtained with re- 
cently developed density-dependent meson-exchange ef- 
fective interaction DD-ME2 [62] , which represents a new 
class of the RMF parametrizations as compared with NLl 
and NL3. However, so far this interaction has not been 
used in the studies of rotating nuclei, thus, it is not em- 
ployed in the current study since its reliability in the 
description of rotational properties is not known. 



III. HYPERDEFORMATION AT HIGH SPIN: 
WHERE TO EXPECT AND ITS GENERAL 
FEATURES 



FIG. 5: (Color online) The same as in Fig. U but for the Te, 
Sn, and Pd isotopes. 



C. The selection of the RMF parametrization. 

The NLl parametrization of the RMF Lagrangian [5^ 
is used in the majority of the calculations in the cur- 
rent manuscript. As follows from previous studies, this 
parametrization provides a good description of the mo- 
ments of inertia of the rotational bands in unpaired 



A. The systematics of crossing spins and transition 
quadrupole moments of the HD bands 

Figs.|4l[5l[6]and[7|display the spins at which the SD and 
HD configurations become yrast (crossing spins) in the 
CRMF calculations. In addition, the calculated transi- 
tion quadrupole moments of these configurations at spin 
values close to the crossing spins are shown. The cal- 
culated HD configurations are near-prolate. One can see 
that the crossing spins are typically lower for proton- 
rich nuclei. Such a feature is seen in most of the isotope 
chains; by going from the /3-stability valley toward the 
proton-drip line, one can lower by approximately 
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FIG. 6: (Color online) The same as in Fig. |4l but for the Cd 
isotopes. The results of the calculations with the NLl (HD 
- open squares, SD - solid squares) and NL3 (HD - solid tri- 
angles up, SD - open triangles down) parametrizations of the 
RMF Lagrangian are presented. Note that the calculations 
with NL3 were performed only for selected nuclei. 



lOh. The minimum of crossing spins is reached at 
N « Z+10 in the Pd, Te and Ru isotope chains (see Figs. 
[Sj3,[n^ and[7^), and the Mo isotope chain (Fig.[71:) shows 
almost no dependence of on mass number. In other 
isotope chains, the minima in crossing spins appear 
in most proton-rich nuclei. Considering that the sensi- 
tivity of modern 7-ray detectors allows to study discrete 
rotational bands only up to « 65ft in medium mass nuclei 
[S^, m, UHl , and that the observation of higher spin states 
will most likely require a new generation of 7-ray track- 
ing detectors such as GRETA or AGATA, these features 
of crossing spins I^.^ represent an important constraint. 

As suggested by the studies of the Jacobi shape tran- 
sition in Ref. [3l|, the coexistence of the SD and HD 
minima at the feeding spins may have an impact on the 
survival of the HD minima because of the decay from the 
HD to SD configurations. If this mechanism is active, 
then only the nuclei in which the HD minimum is lower 
in energy than the SD one at the feeding spin and/or the 
nuclei characterized by the large barrier between the HD 
and SD minima will be the reasonable candidates for a 
search of the HD bands. Figs. [H [5l [6] and [7] show that 
the HD configurations become yrast at lower spin than 
the SD ones only in a specific mass range which depends 
on the isotope chain. This range can be narrow as in the 
case of Te isotopes (Fig. [5^) or wide as in the case of Ce 
isotopes (Fig. 2^). The question of the population of the 
HD bands within the RMF framework definitely deserves 
an additional study, but such a study is beyond the scope 
of the present manuscript. 

Fig. [S] compares the results of the calculations for Cd 
isotopes obtained with the NLl and NL3 parametriza- 



tions of the RMF Lagrangian. One can see that both 
parametrizations predict similar crossing spins I^^P and 
and similar transition quadrupole moments. How- 
ever, in average, the crossing spins calculated with 
NL3 are somewhat lower (by 1 — 2h) than the ones ob- 
tained in the calculations with NLl. 



B. The A ~ 120 region: the analysis of 
experimental data 



TABLE 1: The values of the dynamic moment of inertia Jexp 
of ridge structures measured in the HLHD experiment [2^ . 
Theoretical results obtained in the MM calculations [sj] are 
shown in the last column. 

Nucleus Ji% 



i26Ba 


77 


118 


^23Xe 


71 




i22Xe 


77 


108 


i^iXe 


63 




120 Te 


71 




iisrpe 


111 


97 


i^^Cs 


100 


106 


i^^Cs 


111 




i24xe 


111 


111 


122j 


71 




121j 


77 


102 


i26Xe 


83 


110 



Recent Hyper-Long-HyperDeformed (HLHD) experi- 
ment at the EUROBALL-IV 7-detector array revealed 
some features expected for HD nuclei [13, [32, [s^l ■ Al- 
though no discrete HD rotational bands have been iden- 
tified, rotational patterns in the form of ridge-structures 
in three-dimensional (3D) rotational mapped spectra are 
identified with dynamic moments of inertia J*^^^ rang- 
ing from 71 to 111 MeV~^ in 12 different nuclei selected 
by charged particle- and/or 7-gating (see Table HI. The 
four nuclei, "^Xe, 124^3^ 125^^ ^^^^ 124^6, found with 
moment of inertia J^^) ^ HQ MeV^^ are most likely hy- 
perdeformed [77| while the remaining nuclei with smaller 
values of J^^) g^j-e expected to be superdeformed. The 
width in energy of the observed ridges indicates that 
there are « 6 — 10 transitions in the HD cascades, and 
a fluctuation analysis shows that the number of bands 
in the ridges exceeds 10. The HD ridges are observed in 
the frequency range of about 650 to 800 keV, and their 
dynamic moments of inertia have typical uncertainty of 
10% (e.g. Ill ± 11 McV-i in ^^'^Xe) 

The experimental data show unusual features never be- 
fore seen in the studies of the SD bands. For example, 
the addition of one neutron on going from ^^"^Cs to ^^^Cs 
decreases the experimental J^^) value by 10% (from 
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FIG. 7: (Color online) The same as in Fig. SI but for the Ru, 
Mo and Zr isotopes. 



Ill MeV-i down to 100 MeV'S see Tabled. A similar 
situation is also seen in the SD minimum: the addition of 
one neutron on going from ^^^Xe to ^^^Xe increases the 
experimental J^^) value by - 22% (from 63 MeV-^ to 77 
MeV~^, see Table IJ). It is impossible to find an expla- 
nation for such a big impact of the single particle on the 
properties of nuclei: previous studies in the SD minima 
in different parts of the nuclear chart never showed such 
features. The case of the pair of ^^^Xe and ^^^Xe is even 
more intriguing: a single particle triggers the transition 
from the SD to HD minima (see Table [l|. Considering 
the fact that the ridges corresponding to the SD and HD 
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FIG. 8: (Color online) Calculated kinematic and dynamic mo- 
ments of inertia (top panels) and transition quadrupole mo- 
ments (bottom panels) as a function of rotational frequency 
for the lowest HD solutions in "*Te, ^^^'^^^Cs and ^^^Xe. The 
structure of calculated configurations is indicated at bottom 
panels. Experimental data for dynamic moments of inertia of 
ridge structures are shown in top panels. 




Angular momentum I(fi) 

FIG. 9: (Color online) Energies of the calculated configu- 
rations relative to a smooth liquid drop reference AI[I -\- 1), 
with the inertia parameter A = 0.01. The ND and SD yrast 
lines are shown by dotted and dot-dot-dashed lines, respec- 
tively. Solid and dot-dashed lines are used for the [1,2] and 
[1,21] HD configurations, respectively. Dashed lines represent 
excited HD configurations. 



minima are observed in neighboring nuclei, it is difficult 
to understand why the ridges corresponding to both min- 
ima have not been seen in the same nucleus. 

The calculated kinematic and dynamic moments of in- 
ertia as well as transition quadrupole moments of the low- 
est HD solutions in the candidate HD nuclei are shown 
in Fig. [51 The calculated J^^^ moments of inertia some- 
what underestimate experimental data. The results of 
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the MM calculations for "«Te, i24xe and ^^scs (see Ta- 
ble |T]) are closer to experimental data, but they are ob- 
tained at fixed quadrupole deformation /?2 while other 
deformation parameters f]^, (3q and /Jg are automatically 
readjusted so as to minimize the total free Routhian for 
the vacuum configuration. 
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FIG. 10: Proton (top panel) and neutron (bottom panel) 
single-particle energies (routhians) in the self-consistent ro- 
tating potential as a function of the rotational frequency Q^- 
They are given along the deformation path of the yrast HD 
configuration (the [1,2] conf in Fig. |9]) in ^'^'^Xe and ob- 
tained in the calculations with the NLl parametrization of 
the RMF Lagrangian. Long-dashed, solid, dot-dashed and 
dotted lines indicate (ir = +,r — (vr = +,r — —i), 

(tt = — ,r = +i) and {tt = ~,r — ~i) orbitals, respectively. 
At Q.X ~ 0.0 MeV, the single-particle orbitals are labeled by 
the asymptotic quantum numbers [Nn^ A]Q (Nilsson quantum 
numbers) of the dominant component of the wave function. 
Solid (open) circles indicate the orbitals occupied (emptied) 
in the [1,2] configuration. The dashed box indicates the fre- 
quency range corresponding to the spin-range / = 60 — 85h 
in this configuration. 

In the MM calculations, the kinematic moments of in- 
ertia of the configurations in the HD minimum decrease 



smoothly with the spin, while their dynamic moments of 
inertia are nearly constant (see Figs. 10 and 11 in Ref. 
[sUp. The behaviour of these observables as a function 
of rotational frequency (or spin) is completely different 
in the self-consistent CRMF calculations (see Figs. |8l 
111! and Fig. 1151 below). The kinematic moment of iner- 
tia is either nearly constant or very gradually increases 
with rotational frequency. The dynamic moment of in- 
ertia gradually increases over the calculated frequency 
range showing the features typical to the SD bands in 
the A ~ 190 mass region which are affected by pairing 
d, [S^l : this is despite the fact that pairing is neglected 
in the CRMF calculations. The transition quadrupole 
moment Qt is also increasing with rotational frequency; 
such a feature has not been seen before in the calcula- 
tions without pairing for the SD bands. The microscopic 
origin of these unusual features will be discussed in more 
details in Sect. lmCl 



C. ^^^Xe nucleus 

The results of the CRMF calculations for some HD 
configurations in ^^"^Xe are displayed in Fig. [51 The HD 
minimum becomes lowest in energy at spin 82h, and the 
[1,2] configuration is the yrast HD configuration in the 
spin range of interest. The occupation of the single- 
particle orbitals in this configuration is presented in Fig. 
[TOl The excited HD configurations displayed in Fig. [9] 
are built from this configuration by exciting either one 
proton or one neutron or simultaneously one proton and 
one neutron. The total number of excited HD configu- 
rations shown is 35. It interesting to mention that the 
configuration involving the lowest N — 8 neutron orbital 
(the [1,21] conf. in Fig. [H]) is calculated at low excitation 
energy. 

The calculations reveal a high density of the HD con- 
figurations which will be even higher if the additional 
calculations for the excited configurations would be per- 
formed starting from the low-lying excited HD configura- 
tions, such as the [1,21] configuration. This high density 
is due to two facts: relatively small Z — 54 and iV = 70 
HD shell gaps in the frequency range of interest (see Fig. 
[TO]) and the softness of the potential energy surfaces in the 
HD minimum. Fig.lllb illustrates the latter feature: the 
particle-hole excitations discussed above, characterised 
by low excitation energy, lead to appreciable changes in 
the transition quadrupole moments Qt- It is interesting 
to mention that there are large similarities between the 
single-particle routhians in the vicinity of the Z = 54 
and TV = 70 HD shell gaps obtained in the CRMF cal- 
culations for yrast HD configuration in ^^^Xe (Fig. [TO]) 
and the ones obtained in the Woods-Saxon calculations 
for the HD minimum in ^^^Xe employing the so-called 
universal parametrization of the Woods-Saxon potential 
(see Figs. 8 and 9 in Ref. [3l|). As a consequence, the 
high density of the excited HD states in "'^^''Xe is also ex- 
pected in the MM calculations based on the formalism of 
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FIG. 11: (Color online) Dynamic moments of inertia j''^-' 
(panel (a)) and transition quadrupole moments Qt (panel (b)) 
of the HD configurations in ^^''Xe shown in Fig. [9l They 
are displayed as a function of rotational frequency Qx- The 
regions of band crossings are excluded in these plots. 

The high density of the HD configurations may ques- 
tion our neglect of pairing. This is because there are nu- 
merous possibilities to scatter proton and neutron pairs 
and this process is energetically inexpensive due to the 
high density of the calculated configurations. In order to 
test the impact of pairing on the moments of inertia and 
binding energies, the comparative studies of the vacuum 
HD configuration and its unpaired analog in ^^^Xe and 
of the vacuum SD configuration and its unpaired analog 
in ^^^Dy have been performed within the cranked rel- 
ativistic Hartree-|-Bogoliubov (CRHB) [Hi and CRMF 
approaches. An approximate particle number projection 
by means of the Lipkin-Nogami method is employed in 
the CRHB approach. Note that unpaired analog of the 
vacuum HD configuration in ^^"^Xe (built from the [1,2] 
configuration by the excitation of the proton from the 
7r[770]l/2(r = +i) orbital into the 7r[420]l/2(r = +i) 
orbital, see Fig. I10|) is non-yrast in the spin range of in- 
terest. As follows from this study, in both nuclei the 
pairing has a similar impact on the moments of inertia of 
the configurations under consideration. Taking into ac- 
count that the SD bands in the A ~ 150 mass region are 
well described in the calculations without pairing [1, |3^ , 
it is reasonable to expect that the neglect of pairing is a 
valid approximation for the moments of inertia of the HD 
bands in ^^''Xe. Pairing leads to an additional binding 
of ~ 500 keV in the case of yrast SD band in ^^^Dy; this 
additional binding slightly exceeds 1 MeV in the case of 
the vacuum HD configuration in ^^'*Xe. The dominant 
effects in the quenching of pairing correlations are the 
Coriolis antipairing effect and the quenching due to shell 
gaps: the latter effect being more pronounced in the SD 
bands of the A '-^ 150 mass region because of the larger 



size of the SD shell gaps (see Fig. 4 in Ref. Q). The 
third mechanism of the decrease of pairing is the block- 
ing effect [6^ . Due to this effect the impact of pairing on 
physical observables will be even lower in the HD bands 
of ^^^Xe based on the excitation(s) of one (two) particles 
considered in Fig. [51 Thus, although weak pairing will 
somewhat modify the relative energies of different config- 
urations, in no way will it create an energy gap between 
the vacuum and excited configurations. 
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FIG. 12: The self-consistent neutron density p„{y,z) as a 
function of y- and z- coordinates for the [1,2] configuration 
in ^^''Xe at rotational frequency = 0.75 MeV. Top and 
bottom panels show 2- and 3-dimensional plots of the density 
distribution, respectively. In the top panel, the densities are 
shown in steps of 0.01 fm~^ starting from pn{y,z) — 0.01 
fm~^ 



The calculations suggest that it will be difficult to ob- 
serve discrete HD bands in ^^^Xe since their high density 
will lead to a situation in which the feeding intensity 
will be redistributed among many low-lying bands, thus 
drastically reducing the intensity with which each indi- 
vidual band is populated. On the other hand, the high 
density of the HD bands may favor the observation of 
the rotational patterns in the form of ridge-structures 
in three-dimensional rotational mapped spectra as it has 
been seen in the HLHD experiment [2^ . 

Fig. [8] shows that the HD shapes undergo a centrifu- 
gal stretching that result in an increase of the transi- 
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FIG. 13: (Color online) Kinematic (J<^') and dynamic (j'^') 
moments of inertia as well as transition quadrupole Qt and 
mass hexadecupole Q^o moments of the [1,2] configuration in 
'^^''Xe calculated with different parametrizations of the RMF 
Lagrangian. 



tion quadrupole moments Qt with increasing rotational 
frequency. This process also reveals itself in the mo- 
ments of inertia: the kinematic moments of inertia are 
either nearly constant or slightly increase with increas- 
ing rotational frequency, while the dynamic moments of 
inertia increase continuously and substantially over the 
frequency region of interest. On the contrary, the dy- 
namic moments of inertia of the HD bands are almost 
constant as a function of rotational frequency in the MM 
calculations (see Figs. 10 and 20 in Ref. [311 )j which is 
most likely a consequence of fixed quadrupole deforma- 
tion. The above mentioned features are general ones for 
the HD bands in the A ^ 120 mass region, see Figs. HI 
[TT] and [151 They are in complete contract to the fea- 
tures of the SD bands in unpaired regime, in which the 
Qt, J*-^-* and J'^^-' values (apart from the unpaired band 
crossing regions) decrease with increasing rotational fre- 
quency (see Refs. [4, 37, 44,, ^4£] and references therein). 

Systematic analysis of the yrast/near-yrast HD config- 
urations in the part of the nuclear chart under investi- 
gation shows that the centrifugal stretching is a general 
feature. At the spins, where the HD minimum is lowest in 
energy, it reveals itself (with very few exceptions) by the 
increase of transition quadrupole Qt and mass hexade- 
capole Q^Q moments. Only in a few HD bands, charac- 
terized by the modest transition quadrupole moment, at 
low rotational frequencies these quantities decrease with 
increasing fJ^;. However, even in these bands the Qt and 
Qio values start to increase above specific value of rota- 
tional frequency. Similar features are also seen in the dy- 
namic moments of inertia; with a few exceptions the J^^^ 
values increase in the spin range of interest. The varia- 
tions (both the increases and decreases) in the kinematic 
moments of inertia are rather small (~ 2% of absolute 
value) in the frequency range of interest. 

The basis of the CRMF model is sufficiently large to 
see if there is a tendency for the development of necking. 
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FIG. 14: The weights a^y of different A'^-components in the 
structure of the wave functions of the indicated orbitals. They 
are shown as a function of rotational frequency. For simplic- 
ity, the region of the crossing between the :/[880]l/2~ and 
!/[411]3/2~ orbitals at ~ 0.55 MeV is removed; dotted 
lines are used in panel (b) to connect the weights correspond- 
ing to the //[880]l/2~ orbital before and after crossing. 



Fig. [T^l shows some indications of the necking and the 
clusterization of the density into two fragments in the 
[1,2] configuration of ^^^Xe, but this effect is not very 
pronounced in this nucleus. 

The kinematic and dynamic moments of inertia as 
well as the transition quadrupole and mass hexadecapole 
moments of the [1,2] configuration in ^^*Xe are shown 
for different parametrizations of the RMF Lagrangian 
in Fig. 1131 The gradual increase of all physical observ- 
ables is due to centrifugal stretching. The NLZ (NLSH) 
parametrizations provide the largest (smallest) values of 
the above mentioned physical observables, while the re- 
sults obtained with NLl and NL3 are in between those 
results. Similar relations between the results obtained 
with these parametrizations also exist in other regions 
of nuclear chart studied so far in the CRMF or CRHB 
frameworks, namely, in the ^ 60 [T^I, A ^ 150 [3§] 
and A ^ 190 m&k regions of superdeformation and in 
the A ~ 250 [54] region of normal deformation. The 
NLl and NL3 parametrizations, which have been exten- 
sively used in the previous studies of rotating systems 
and superdeformation [4^ . give the values of physical 
observables of interest which differ only by few %. It is 
known that the NLSH parametrization somewhat under- 
estimates the experimental moments of inertia [s^, [t^] . 
The NLZ parametrization has not been used in the pre- 
vious studies of rotating systems, so it is unknown how 
well it describes such systems. 
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FIG. 15: (Color online) Dynamic moments of inertia j'^' of selected configurations in ^^''Xe and neighbouring nuclei. Dynamic 
moment of inertia of the [1,2] configuration A in ^^''Xe is shown by a thick solid line in each panel. The J^'^-' values of the 
configurations in the nucleus indicated on the panel are displayed by the lines of other types. These configurations differ from 
the [1,2] configuration A in ^'^'^Xe in the occupation of the orbitals shown in the panels. Vertical dashed lines indicate the 
frequency range corresponding to the spin range 7 = 60 — 85fi in the [1,2] configuration of ^^*Xe. 



D. Single-particle properties at hyperdeformation: 
an example of neighbourhood of ^^^Xe. 

The role of the single-particle degrees of freedom at 
hyperdeformation was mainly overlooked in the previous 
studies. It has been studied to some extent only within 
the MM method in Refs. [HIIH. However, the studies 
of Ref. 31] suggest that the ^^''Xe nucleus is very rigid 
in the HD minimum: the dynamic moments of inertia of 
different HD bands differ by no more than 2%, and their 
changes as a function of spin are very small (see Fig. 10 
in Rcf. [31]). Similar results were obtained for HD bands 
in i46Gd and ^^aoy in Ref. 

On the contrary, the CRMF calculations for the dy- 
namic moment of inertia of the yrast and excited HD con- 
figurations in ^^"^Xe show much larger spread and much 
larger variations as a function of rotational frequency, 
see Fig. [TTb . In addition, large variations in the calcu- 
lated transition quadrupole moments Qt of these con- 
figurations are clearly seen in Fig. [TTb . This suggests 
that the HD minimum is relatively soft and that the in- 
dividual properties of the single-particle orbitals play an 
important role in the definition of the properties of the 
HD bands. One of our goals is to investigate the impact 
of the particle in a specific single-particle orbital on the 
properties of the HD bands and to study whether the 
methods of configuration assignment based on the rela- 



tive properties of different bands are also applicable at 
HD. 



1. The structure of the wave function 

The structure of the wave function at HD is analysed 
on the example of a few single-particle orbitals of the [1,2] 
configuration in ^^^Xe (Fig. [TT)l . The evolution of these 
orbitals in energy with rotational frequency is displayed 
in Fig. [TOl The wave function 5* is expanded into the 
basis states by 

^ = Y,('NjNa> (9) 

N,a 

where N and a represent the principal quantum number 
and the set of additional quantum numbers specifying 
the basis state, respectively. We specify the weight aj^ 
of the basis states belonging to the specific value of N in 
the structure of the wave function as 

E (10) 

A'' — fixed, a 

with the condition aj^ — 1 following from the or- 
thonormalization of the wave function of the single- 
particle orbital. 
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FIG. 16: (Color online) Effective alignments i^ff extracted from the calculated configurations for the orbitals active in the 
vicinity of the Z = 54/55 and N = 70 HD shell gaps (see Fig. llUp . The calculated configurations are the [1,2] conf. in ^^*Xe and 
the configurations in neighboring nuclei (shown in Fig. I15|) obtained by adding or removing a single particle (proton or neutron). 
The effective alignment between configurations X and Y is indicated as "X/Y". The configuration X in the lighter nucleus is 
taken as a reference, so the effective alignment measures the effect of the additional particle. The compared configurations 
differ in the occupation of the orbitals shown in the panels. Note that the vertical scale of different panels is different. Vertical 
dashed lines indicate the frequency range corresponding to the spin range / = 60 — 85?i in the conf. A of ^^^Xe. 



Hyperdeformation leads to a considerable fragmenta- 
tion of the wave function over N, which is much larger 
than in the case of SD. In the regions away from the band 
crossing the weight aj^ of the dominant A^-component of 
the wave function does not exceed 0.8 while the weight 
of second largest component is typically around 0.2 (Fig. 
I14p . Very strong fragmentation of the wave function is 
seen in the case of the z^[761]3/2+ orbital: before the 
band crossing the weights of the N = 7 and = 5 com- 
ponents of the wave function are approximately 0.6 and 
0.3, respectively. Even stronger fragmentation is seen in 
the region of the band crossing of the [761] 3/2+ and 
:/[301]3/2+ orbitals at fl^ ~ 0.7 MeV (Figs. [TO]) where 
they strongly interact and gradually exchange their char- 
acter (Figs. [Hk and c). Similar fragmentation is also 
seen for the 7r[770]l/2+ orbital (Fig. [Ti|) which interacts 
strongly with the tt [532] 5/2+ orbital in the band crossing 
region at - 0.8 MeV (Fig. [lOl). 



2. The methods of configuration assignment 

The HD bands in nuclei neighboring to ^^^''Xe, which 
differ by either one proton or one neutron from the [1,2] 
configuration in ^^^^Xe, and their relative properties with 
respect of the [1,2] configuration in ^^^Xe are studied in 
order to investigate the applicability of different methods 
of configuration assignment at HD. 

The dynamic moments of inertia for the four HD bands 
in each of these nuclei are compared with the one of the 
[1,2] configuration in ^^^^Xe in Fig. [151 The difference 
between the dynamic moments of inertia of the configu- 
rations in nuclei with masses A and A ± 1 is due to the 
impact of the particle in the specific single-particle or- 
bital by which two compared configurations differ. The 
results of the calculations question conventional wisdom 
[stI ] that the largest impact on the dynamic moment of 
inertia is coming from the particles in the intruder or- 
bitals. Indeed, the impact of the neutron in the hyperin- 
truder i/[880]l/2~ orbital on the dynamic moments of in- 
ertia (Fig. [TsH) is comparable to the one of non-intruder 
[642] 5/2+ orbital or even smaller by a factor of ~ 2 than 
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FIG. 17: (Color online) Relative transition quadrupole moments AQt = Qt{A -\- 1) — Qt{A) [A is the mass of the nucleus] 
extracted from the calculated configurations in indicated nuclei. The compared configurations are shown as "X/Y": the 
configuration X in the lighter nucleus is tal<en as a reference, so the AQt measures the effect of the additional particle placed in 
the orbitals shown in the panels. Vertical dashed lines indicate the frequency range corresponding to the spin range / = 60 — 85/1 
in the [1,2] configuration of ^■^'*Xe. 



tlie impact due to tlie neutron in non-intruder j/[532]3/2+ 
orbital (Fig. [TSbV A similar situation is also seen for 
protons, where, for example, the impact of the proton in 
the hyperintruder 7r[770]l/2+ orbital is smaller than its 
impact in the non- intruder tt [420] 1/2^ orbital. This sug- 
gests that not only angular momentum, carried by the 
particle in specific single-particle orbital, but also po- 
larization effects it induces into time-even and time-odd 
mean fields (sij are important when considering relative 
properties of two configurations. Based on this example, 
one can conclude that the configuration assignment of 
the HD bands, based only on the relative properties of 
the dynamic moments of inertia of two compared bands, 
is unreliable. 

The configuration assignments at SD have been mostly 
based on the effective alignment approach (see Refs. 
[H, m, IZ^I and references therein) . The success of this 
method is due to the fact that it was possible to separate 
intruder and non-intruder orbitals since the former show 
pronounced dependence of the effective alignments ieff 
on the rotational frequency (see, for example. Figs. 2, 3, 
5, 6 and 8 in Ref. [23]). On the contrary, the effective 
alignments of non-intruder orbitals are typically constant 
as a function of rotational frequency. It also follows from 



the studies in the A ^ 140 — 150 region of superdefor- 
mation that the change of effective alignment by « l/i 
within the observed frequency range allows to identify 
aligning intruder orbitals with a high level of confidence. 

A configuration assignment based on the effective 
alignments depends on how accurately these alignments 
can be predicted. For example, the application of the ef- 
fective alignment approach in the A ^ 140—150 region of 
superdeformation requires an accuracy in the prediction 
of if,ff on the level of Q.ih and ~ 0.5/i for nonintruder 
and intruder orbitals, respectively [3, [H, [1^. In the 
highly deformed and SD bands from the A ~ 60 — 80 mass 
region, these requirements for accuracy are somewhat re- 
laxed (HjIt^]. We expect that in the A ~ 125 mass region 
of HD, the effective alignments should be predicted with 
a precision similar to that in the A ^ 140 — 150 region 
for a reliable configuration assignment. 

Our analysis shows that a reliable configuration as- 
signment for the HD bands based solely on the effective 
alignment approach will be problematic (at least in the 
A ~ 125 mass region) because of several reasons. First, 
the hyperintruder orbitals do not show appreciable varia- 
tions of ief f with rotational frequency. Fig. ll6l shows that 
the effective alignments of the hyperintruder orbitals such 
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FIG. 18: (Color online) Proton and neutron single-particle energies in ^"^Cd, as a function of charge quadrupole moment 
Q, obtained in the axially symmetric RMF calculations. Solid and dashed lines denote positive and negative parity orbitals, 
respectively. The Fermi energy Ep is shown by dotted line. The single-particle orbitals are labeled by the Nilsson quantum 
numbers. Large shell gaps are indicated. 



as 7r[770]l/2"'" and i^[880]l/2~ show little variations with 
rotational frequency (see Fig. 116b . dV On the contrary, 
the effective alignments of the :^[532]3/2~'" and i^[530]l/2^ 
orbitals show much larger variations reaching l.bh in the 
spin range / = 60 — 85ft in the case of the latter orbital 
(see Fig. 116b). However, the variations of ig// as a func- 
tion of rotational frequency are small for the majority 
of the orbitals in the spin range of interest. Thus, con- 
trary to the case of SD, it will be more difficult to distin- 
guish between hyperintruder, intruder and non-intruder 
orbitals based on the variations of Zg// with rotational 
frequency. This situation will become even more compli- 
cated if the suggestion of Ref. [3l| that the spin range 
over which the HD bands are expected to be observed 
{2Ah at the most; this is shorter than in the case of SD) 
is true. These two features (small variations of ieff and 
expected spin (frequency) range of the HD bands) will 
lead to a situation where the ieff values for many or- 
bitals will look alike within the typical 'error bars' of 
the description of ieff by theoretical models, so that it 
will be difficult to distinguish between them within the 
framework of the effective alignment approach. 



Similar to the case of SD [i^ |4l| , additional informa- 
tion on how the single particle affects the properties of 
the HD bands can be extracted from the relative transi- 
tion quadrupole moments AQt . Fig. [T7] shows that the 
hyperintruder 7r[770]l/2+ and i^[880]l/2- orbitals with 
AQt « 2 eb and AQt « 1.25 eb have the largest im- 
pact on the transition quadrupole moments among the 
studied proton and neutron orbitals. One has to keep in 
mind that the addition of a proton changes the proton 
number by one. This change contributes approximately 
0.5 eb in relative transition quadrupole moment AQ* of 
the proton orbitals. This effect is not present in the AQt 
values of the neutron orbitals. 

The AQt values were used only as a complimentary 
tool of the configuration assignment at SD. This is be- 
cause of the difficulty to measure them in experiment 
72, 73] and the fact that they show little variation as a 
function of rotational frequency, thus providing less in- 
formation than ieff - The same features are also valid at 
HD; see Fig. [T7] for the variations of the AQt values. In 
addition, some single-particle orbitals such as 7r[422]3/2^ 
and 7r[303]7/2^ (Fig. [T7b) show very similar AQj values. 
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FIG. 19: The energy gaps between the last occupied and first 
unoccupied single-particle orbitals shown as a function of neu- 
tron number for different isotope chains. They are extracted 
from the routhian diagrams of the lowest HD configurations 
at the spin values where these configurations become yrast. 
The bars are used for the energy gaps in the Cd isotopes. 



This will not allow to make a unique configuration as- 
signment even if the experimental AQt values for these 
orbitals are available. On the other hand, their Zg// val- 
ues differ by ~ IH (Fig. [T6h). and this fact can be used 
in the configuration assignment. 

However, the fact that in general the effective align- 
ment approach fails to provide a unique configuration 
assignment at HD increases the role of the method of 
configuration assignment based on relative transition 
quadrupole moments. Our analysis shows that only si- 
multaneous application of these two methods by compar- 
ing experimental and theoretical (ieff,A.Qt) values will 
lead to a reliable configuration assignment at HD. 

Let us illustrate this on the hypothetical example of 
two "experimental" bands; one in ^^'^I and another in 
^^''Xe. In this example, the [1,2] configuration is assigned 
to the band in ^^''Xe. Let us assume that the effective 
alignments in the ^^'^I/^^^'Xe pair of the bands increase 
from i.Oh to 4:.25h in the frequency range 0.62-0.87 MeV 



under selected spins of these bands. Under these con- 
ditions, the "experimental" bands differ in the occupa- 
tion of the 7r[770]l/2+ orbital (Fig. [16ji). However, it 
is reasonable to expect that the spins of "experimental" 
bands will not be fixed, so these changes in effective align- 
ment should be from (4.0 + n)h to (4.25 -I- n)h, where 
n — 0, ±1,±2,.... Assuming that the accuracy of the 
description of effective alignments in theoretical calcula- 
tions is around 0.4?i, one can conclude that for n = —3 
the "experimental" bands can also differ in the occupa- 
tion of either the 7r[532]5/2~ or 7r[651]3/2+ orbitals (Fig. 

t). In a similar way to the A ^ 150 region of SD 
[39| , the systematic studies of the pairs of the bands 
which differ by one proton may narrow the choice of the 
orbitals involved. On the other hand, the AQt values for 
these orbitals are drastically different; AQt ~ 2.0 eb for 
the 7r[770] 1/2+ orbital, AQt « 1.4 eb for 7r[651]3/2+, and 
AQt « 0.7 eb for 7r[532]5/2- (see Fig. HT]). So, if both 
quantities, ieff and AQt, are measured simultaneously, a 
unique configuration assignment for "experimental" band 
in ^^"^I will be possible. 

The band crossing features of the HD bands provide 
an additional tool of configuration assignment which can 
be used more frequently than in the case of the SD bands 
because of strong mixing between the different A^-shells 
at HD. The large peaks in J^^^ of the and vB configu- 
rations in ^^^Xe (Fig. fTsH) are due to the band crossings 
with a strong interaction. These crossings are also visible 
in the effective alignments ieff (Fig. \W\i) and relative 
transition quadrupole moments AQt (Fig. I17H ). They 
originate from the crossing of the same signatures of the 
;/[301]3/2 and i/[761]3/2 orbitals, where i^A and i/B have 
signatures r — -\-i and r — —i, respectively. The former 
orbital is occupied before band crossing, the latter after 
band crossing. An unusual feature of these band cross- 
ings is the fact that they originate from the interaction of 
the orbitals, the dominant A^-components of which differ 
by AA^ = 4. At SD, the crossings between the orbitals 
dominated by different A^-shells have been characterized 
by a weak interaction leading to a sharp jump in J*-^-* 
0, III, [ill- The observed unpaired SD band crossings 
with strong interaction are between the orbitals with the 
same dominant A'^-shells and they were observed in the 
nuclei around "^Gd 



E. 



General observations: the density of the HD 
bands and the necking degree of freedom 



As discussed in Sect. IIII CI on the example of ^^''Xe, 
the high density of the HD bands is one of the major 
obstacles for the observation of discrete HD bands. It 
will lead to a situation where the feeding intensity will 
be redistributed among many low-lying HD bands, thus, 
drastically reducing the intensity with which each indi- 
vidual band is populated. As a consequence, the feeding 
intensity of an individual HD band will drop below the 
observational limit of experimental facility; this fact has 
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to be taken into account when planning future experi- 
ments for a search of discrete HD bands. 

Two factors contribute to the high density of the HD 
bands, namely, relatively small proton and neutron HD 
shell gaps in the frequency range of interest and the soft- 
ness of the potential energy surfaces in the HD minimum 
(see Sect. IIII CI) . Systematic mapping of the density of 
the HD states as a function of the proton and neutron 
numbers is too costly in the computational sense because 
it involves the calculation of the lowest in energy particle- 
hole excitations. Thus, we decided to look at the problem 
of the density of the HD states in a somewhat simplis- 
tic way by considering the proton and neutron energy 
gaps between the last occupied and the first unoccupied 
states in the yrast HD configurations; the small size of 
these gaps will most likely point to the high density of 
the HD bands. 

The analysis of the Nilsson diagrams in Fig. [15] already 
reveals some HD gaps in the single-particle spectra. At 
the values of Qo ^ 17 — 20 eb typical for the HD con- 
figurations in Cd isotopes (Fig. [BJs) , there are very large 
proton Z = 48 and neutron = 48 HD shell gaps and 
smaller neutron gaps at = 58 and 60. In general, this 
figure suggests that the hyperdeformation will be more 
favoured in the nuclei with a similar number of protons 
and neutrons because the proton and neutron shell effects 
for the HD shapes will act coherently; this trend has al- 
ready been seen in the crossings spins for different 
isotope chains in Sect. IIII Al 

The size of these gaps and their presence will be altered 
(especially, for medium and small size energy gaps) when 
the rotation and the self-consistent readjustment of the 
neutron and proton densities with the change of particle 
number are taken into account. Indeed, this is seen in 
Fig. [Tni which shows the energy gaps between the last 
occupied and first unoccupied single-particle orbitals as 
a function of the neutron number for different isotope 
chains. The largest proton gap at Z = 48 is seen in 
Cd isotopes; its size is around 1.5 MeV in proton-rich 
nuclei and it increases up to 3 MeV with the increase of 
neutron number. In other isotope chains, the size of the 
proton energy gap is smaller than in Cd isotopes and it 
fluctuates around 1 MeV. For the majority of the nuclei, 
the size of the neutron energy gap fluctuates around 1 
MeV. However, its size increases up to 1.5 MeV in some 
nuclei and in ^^Cd it reaches 2 MeV (see Fig. [12] for 
details) . 

Taking into account that the proton and neutron HD 
shell gaps in ^■^'*Xe are around 1 MeV (Fig. [TU)) and con- 
sidering the results for the density of the HD states in this 
nucleus as a reference fSect- lHTCj) . one can conclude that 
the analysis of the energy gaps suggests that in most of 
the nuclei the density of the HD bands will be high. For 
these nuclei, the observation of discrete HD bands using 
existing facilities is most likely not possible. The only 
exceptions are Cd nuclei and a few nuclei in which the 
size of at least one gap reaches 1.5 MeV (see Fig. [TO] for 
details). For example, in Cd nuclei the large size of the 



Z = A8 HD shell gap (especially, for nuclei in the valley of 
the /^-stability) will make proton particle-hole excitations 
energetically expensive. As a consequence, the density of 
the HD bands has to be lower in Cd isotopes as compared 
with the one in other isotopes. 

One has to remember that the high density of the 
HD bands is not necessarily a negative factor. It favors 
the observation of the rotational patterns in the form of 
ridge-structures in three-dimensional rotational mapped 
spectra as it has been seen in the HLHD experiment for a 
few nuclei [22| . The observation of ridge-structures as a 
function of proton and neutron number, which seems to 
be feasible with existing experimental facilities such as 
GAMMASPHERE, will provide invaluable information 
about HD at high spin. 

The importance of the necking degree of freedom for 
the high-spin HD states has been studied in the MM ap- 
proach in Refs. [1, [IB]- However, this degree of freedom 
has not been investigated in detail at high spin in self- 
consistent approaches so far. In order to flU this gap 
in our knowledge, the systematics of the self-consistent 
proton density distributions in the HD states obtained in 
the CRMF calculations are shown in Fig. [20] One can see 
that in some nuclei such as ^^''Te, ^^°Xe, ^"^^Ba the neck- 
ing degree of freedom plays an important role, while oth- 
ers (for example, ^'^°Mo and ^'^^Ce) show no necking. The 
neck is typically less pronounced in the HD states of the 
lighter nuclei because of their smaller deformation (see 
also Fig. 5 in Ref. ^2^). It becomes even more important 
in extremely deformed structures which according to the 
language of Ref. i5] can be described as megadeformed. 
Fig. [21] shows an example of density distribution for the 
megadeformed state in ^°^Pd, which becomes yrast at 
/ ~ 85S in the CRMF calculations. The neck is more 
pronounced in the proton subsystem than in the neutron 
one both in the HD and megadeformed structures due 
to the Coulomb repulsion of the segments. This is illus- 
trated in Fig.[5T] Our self-consistent calculations indicate 
that the shell structure is also playing a role in a forma- 
tion of neck. For example, the neck is visible in ^^^Ba 
but is not seen in ^^^Ba fFig.l^O]). This is contrary to the 
fact that the calculated transition quadrupole moments 
of the HD states in these nuclei (Fig. [4ji) and their den- 
sity elongations (Fig. [^Ojl are comparable. These results 
indicate that, in general, the necking degree of freedom is 
important in the HD states and that it should be treated 
within the self-consistent approach which, in particular, 
allows different necking for the proton and neutron sub- 
systems. 



IV. "^I NUCLEUS: A CANDIDATE FOR A 
DOUBLY MAGIC EXTREMELY SD BAND. 

The results of the CRMF calculations for the conflgu- 
rations forming the yrast line or located close to it in en- 
ergy are shown in Fig. [22 According to the calculations, 
normal- and highly-deformed bands, many of which show 



18 



116 



Ba 




— 1 1 1 1 r 

132| 



Ba 




— 1 1 1 r- 



106- 



Te 




124- 



Te 




124, 



Ce 




— 1 1 1 r- 



110 



Xe 




136, 



Ce 




i30Xe 




98 



Cd 





104, 



Sn 





92 



Ru 





92 



Pd 




"Pd 




1 1 1 1 1 


1 1 1 1 1 

90Zr 


1 1 1 1 1 

88 Mo 


1 1 1 1 1 

ioOmo 






1 1 1 1 1 




1 1 1 1 1 


1 1 1 1 1 



-10 



10 



-10 



5 10 -10 -5 
Symmetry axis z (fm) 



10 



-10 



10 



FIG. 20: The self-consistent proton density pp{y,z) as a function of y- and z- coordinates for the HD configurations. They 
are displayed at spin values at which these configurations become yrast. For each isotope chain, the densities in two nuclei 
(typically, most proton- and neutron rich ones included in calculations) are shown. The densities are displayed in steps of 0.01 
fm""^ starting from pp{y,z) = 0.01 fm~^. 



the high triaxiality that is indicative of approaching band 
termination ^55i |. dominate the yrast hne up to / w 64?i. 
At higher spin, more deformed structures become yrast. 
The configuration A has the structure ttQ^i/Q^ and is yrast 
in the spin range / = 64 — 73h: no hyperintruder N = 7 
orbitals are involved in its structure. In this spin range 
it is characterized by the transition quadrupole moment 
Qt ~ 15.7 eb and by the 7-deformation of ~ 1°. The 



normahzed transition quadrupole moment in this system 
is Q"°™ — 11.7 eb, thus, this band is approximately 
35% more deformed than the SD band in ^^^Dy. As 
a consequence, in terms of deformation, this band can 
be characterized as an extremely superdeformed (ESD) 
band which is only slightly less deformed than the HD 
bands. 

In addition, the configuration A is well separated from 
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FIG. 21: The same as in Fig. 1121 but for yrast megadeformed 
state in ^°^Pd at rotational frequency 0,^ = 0.95 MeV. Two 
top panels show 2-dimensionaI plots of the proton and neutron 
density distribution. 



the excited SD/HD configurations below / ^ 73h (see 
Fig. [22]) . This is due to the presence of the large Z — 53 
and = 58 ESD shell gaps in the single-particle spec- 
tra (see Fig. [221) ■ In this configuration, all single-particle 
states below the Z = 53 and iV = 58 ESD shell gaps are 
occupied by protons and neutrons, respectively. Thus, 
this ESD band is a doubly-magic one. This band appears 
as doubly- mag ic also in the calculations with widely 
used NL3 [58| | and NLZ [6l| parametrizations of the 



FIG. 22: (Color online) Energies of the calculated configu- 
rations relative to a smooth liquid drop reference AI{I + 1), 
with the inertia parameter A = 0.01. Normal deformed (ND), 
SD and HD configurations are shown by dotted, dot-dashed 
and solid lines, respectively. Configuration A is shown by 
long-dashed line. 



TABLE II: The size of the Z=53 and N=58 ESD shell gaps 
[in MeV] obtained with different parametrizations of the RMF 
Lagrangian for the configuration A in I at spin / — 60ft 
(rotational frequency ^Ix ~ 0.96 MeV). 
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RMF Lagrangian, see Table HIl Extensive calculations 
with the NL3 parametrization (similar to the ones pre- 
sented in Fig. [52) show that this band become yrast at 
/ - 62h. The Z = 53 ESD shell gap is smaller than 1 
MeV only in the NLSH [6^ parametrization of the RMF 
Lagrangian (see Table [IT]). However, it is known that 
the single-particle energies are not well described in this 
parametrization [53|. One should note, however, that 
the size of the ESD gaps in the configuration A of ^^^I is 
somewhat smaller than the one for the yrast SD band in 
152 j3y (compare Fig. [211 in the present manuscript with 
Fig. 3 in Ref. see also Figs. 4, 11, 12 in Ref. 
obtained with different parametrizations of the RMF La- 
grangian and relevant for ^^^Tb). 

The dynamic moments of inertia of the configuration 
A in ^"'^^I and the configurations in neighboring nuclei 
are shown in Fig. 1241 The increase of J(2) at ~ 1.2 
MeV is in part due to unpaired band crossing caused by 
the interaction of the occupied i/[413]7/2" and unoccu- 
pied z/[651]3/2" orbitals (Fig. 155)) . A centrifugal stretch- 
ing may also contribute to this increase of j'-^K The 
effect of the occupation of a single proton (neutron) in- 
truder orbital on the properties of the ESD bands is much 
more pronounced than that in the HD bands of the nuclei 
around ^^"'Xe (see Sect. IHIDj) : the changes induced into 
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FIG. 23: The same as in Fig. 1101 but for the configuration 
A in ^^^I. Solid (open) circles indicate the orbitals occupied 
(emptied). The dashed box indicates the frequency range 
corresponding to the spin range / = 50 — 75h in this configu- 
ration. 



dynamic moment of inertia reach at least 10% of its abso- 
lute value for the 7r[660] 1/2+ (Fig.[Ml:), 7r[4-f6]l/2+ (Fig. 
[Hi), i/[651]3/2+ (Fig.[2li) and i^[651]3/2- (Fig.[2li) or- 
bitals. In a similar way, the effective alignments of these 
orbitals as well as of the 7r[541]l/2+ orbital show appre- 
ciable variations as a function of rotational frequency (see 
Fig. [25]), reaching at least Ih in the spin range of interest. 
This suggests that the configuration assignment based on 
the effective alignment method will be more reliable in 
the case of BSD bands as compared with the HD bands in 
the nuclei around ^^''Xe (see Sect. IIIIDI for a discussion 
of these methods). Relative properties of the dynamic 
moments of inertia of two compared bands will also play 
a complimentary role in the configuration assignment. 



V. CONCLUSIONS 

For the first time, the hyperdeformation at high spin 
has been studied in a systematic way within the frame- 
work of a fully self-consistent theory: the cranking rel- 
ativistic mean field theory. The study covers even-even 
nuclei in the Z — AO — 58 part of nuclear chart. The main 
results can be summarized as follows: 

• The crossing spins , at which the HD configu- 
rations become yrast, are lower for proton-rich nu- 
clei. This is a feature seen in the most of studied 
isotope chains; by going from the /3-stability valley 
towards the proton-drip line one can lower by 
approximately lOh. 

• The density of the HD bands in the spin range 
where they are yrast or close to yrast is high in 
the majority of the cases. For such densities, the 
feeding intensity of an individual HD band will 
most likely drop below the observational limit of 
modern experimental facilities. This fact has to be 
taken into account when planning the experiments 
for a search of discrete HD bands. Our calcula- 
tions indicate Cd isotopes and few other nuclei with 
large shell gaps (see Sect. IHI El for details) as the 
best candidates for a search of discrete HD bands. 
An alternative candidate is the doubly magic ex- 
tremely superdeformed band in ^^^I, the deforma- 
tion of which is only slightly lower than that of the 
HD bands, and which may be observed with exist- 
ing experimental facilities. 

• The high density of the HD bands will most likely 
favor the observation of the rotational patterns in 
the form of ridge-structures in three-dimensional 
rotational mapped spectra. The study of these pat- 
terns as a function of proton and neutron numbers, 
which seems to be possible with existing facilities, 
will provide a valuable information about hyperde- 
formation at high spin. 

• With a very few exceptions, the HD shapes un- 
dergo a centrifugal stretching that results in an in- 
crease of the values of the transition quadrupole 
Qt and mass hexadecapole Q40 moments as well as 
the dynamic moments of inertia J*-^) with increas- 
ing rotational frequency. The kinematic moments 
of inertia J^^^ show very small variations in the fre- 
quency range of interest. These are general features 
of the HD bands which distinguish them from the 
normal- and superdeformed bands. Such features 
have not been seen before in the calculations with- 
out pairing. In unpaired regime, the Qt, J*-^-* and 
j'^' values decrease with rotational frequency in 
the SD configurations; the only exceptions are the 
regions of unpaired bands crossings. 

• The individual properties of the single-particle or- 
bitals are not lost at HD. In the future, they will 
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FIG. 24; (Color online) The same as in Fig. [151 but for dynamic moments of inertia J^^' of the configurations used in Fig. [25] 
below. Dynamic moments of inertia of the configuration A in ^^^I are shown by solid line in each panel. Vertical dashed lines 
indicate the frequency range corresponding to the spin range / = 50 — 75h in the configuration A of ^^^I. The -k[4 + 6]1/2 label 
in panel (a) indicates the orbital with strong mixing of the = 4 and A'' = 6 shells: this mixing predominantly emerges from 
the interaction of the vr [420] 1/2 and tt [660] 1/2 states. 



allow the assignment of the configurations to the 
HD bands using the relative properties of differ- 
ent bands. Such methods of configuration assign- 
ment were originally developed for superdeforma- 
tion. In contrast to the case of SD, our analysis in 
the A ~ 125 mass region shows that only simulta- 
neous application of the methods based on effective 
alignments and relative transition quadrupole mo- 
ments by comparing experimental and theoretical 
(ie//,A(5t) values will lead to a reliable configu- 
ration assignment for the HD bands. Moreover, 
additional information on the structure of the HD 
bands will be obtained from the band crossing fea- 
tures; the cases of strong interaction of the bands 
in unpaired regime at HD will be more common as 
compared with the situation at SD. 

The physics of hyperdeformation at high spin is also 
defined by the fission barriers; the competition with fis- 
sion certainly makes the population of the HD states dif- 



ficult. It is an important issue, which, however, goes be- 
yond the scope of the current manuscript. It is likely that 
the fission barriers are small or non-existent at the spins 
around 80 — 90h in some of the studied nuclei; the obser- 
vation of the HD bands then will not be possible in these 
systems. This problem definitely deserves a deeper atten- 
tion; the study of the fission barriers at high spin typical 
for HD within the framework of the cranked relativistic 
Hartrec-Bogoliubov theory is in its initial stage and the 
results will be presented in a forthcoming manuscript. 
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